Time-Reversal-Breaking and d-Wave Superfluidity 
of Ultracold Dipolar Fermions in Optical Lattices 
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We describe possible superfluid phases of ultracold dipolar fermions in optical lattices for two- 
dimensional systems. Considering the many-body screening of dipolar interactions at larger filling 
factors, we show that several superfluid phases with distinct pairing symmetries naturally emerge in 
the singlet channel: local s-wave {si), extended s-wave (se), d-wave (d) or time-reversal-symmetry 
breaking {si + se ± id)-wave. The temperature versus filling factor phase diagram indicates that 
d-wave is favored near half-filling, that {si + se)-wave is favored near zero or full filling, and that 
time-reversal-breaking {si + se ± id)-wave is favored in between. When a harmonic trap is included 
a sequence of phases can exist in the cloud depending on the filling factor at the center of the trap. 
Most notably in the region where the {si + se± id)-wave superfiuid exists, spontaneous currents are 
generated, and may be detected using velocity sensitive Bragg spectroscopy. 

PACS numbers: 03.75. Hh, 03.75. Kk, 03.75. Ss, 67.85.-d 
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Ultracold heteronuclear molecules are very interest- 
ing quantum systems to study because they possess 
electric dipolc moments. This internal degree of free- 
dom adds richness to the nature of interactions between 
molecules in comparison to interactions between atoms 
in purely atomic systems. Dipolar molecules can be ei- 
ther fermionic or bosonic in nature depending on their 
constituent atoms, and dipolar interactions allow for the 
emergence of quantum phases which may be extremely 
difhcult to be realized in condensed matter. Recently, 
ultracold dipolar molecules were produced optically [1| 
followed by their production from Bose-Fermi mixtures 
of ultracold atoms first in the vicinity of Feshbach reso- 
nances and later brought into their roto-vibrational 
ground state The production of these heteronuclear 
molecules in harmonic traps has paved the way for stud- 
ies of the quantum phases of interacting dipolar bosonic 
or fermionic molecules. In the case of trapped clouds, 
a few quantum phases have been proposed for dipo- 
lar bosons including ferroelectric superfluids Q, Wigner 
crystals while for dipolar fermions phases such as 
ferroelectric Q or ferro-nematic Fermi liquids and 
Berezinskii-Kosterlitz-Thouless or j-triplet super- 
fluids have been suggested. In the case of optical lat- 
tices, additional phases such as supersolids 0] or micro- 
emulsions [lo| have been proposed for dipolar bosons, 
while studies for dipolar fermions in optical lattices are 
just beginning. 

A natural next step for experiments is the loading of 
dipolar (heteronuclear) molecules in optical lattices. In 
anticipation of these experiments, we discuss here the 
quantum phases of dipolar fermions in harmonically con- 
fined optical lattices, paying particular attention to the 
emergence of superfluid phases that break time reversal 
symmetry spontaneously, as there are no confirmed ana- 
logues in condensed matter physics (Tlj . By including the 
effects of screening at intermediate and high filling fac- 



tors, we argue that quantum phases of dipolar fermions in 
harmonically confined optical lattices (two-dimensional 
square geometry) can be approximately described by an 
extended Hubbard model, with local on-site and a few 
neighbor interactions. For attractive local and nearest 
neighbor interactions, we derive the phase diagram and 
establish all accessible phases in the singlet channel, the 
most important of which corresponds to c?-wave super- 
fluidity and a superfluid phase that breaks time reversal 
symmetry spontaneously and involves a superposition of 
s-wave and d-wave components of the order parameter. 
In the broken time-reversal symmetry phases, sponta- 
neous currents flow and can be detected using velocity 
sensitive Bragg spectroscopy l3, . 

The bare Hamiltonian for dipolar fermions in optical 
lattices for a two-dimensional system (xy-plane) is 



Hh 



where (...) indicates nearest neighbors, and — c-^Cii^ 
is the local particle number operator. The on-site interac- 
tion U = Us + Udip contains two contributions. The first 
is from s-wave scattering Us = {Anh^ as / m) /(ir|'u;(r)|^ 
and the second is from the on-site dipole-dipole inter- 
action Udip = J dk.Wp(k)Vdip(k), where Us is the s- 
wave scattering length, w{r) is the Wannier function, 
wp^k) and Vdip(k) are the Fourier transforms of |'u;(r)p 
and of the dipole-dipole potential Vij, respectively. The 
long-range part of the dipole-dipole interaction is Vij « 
[e; • ej - 3(ej • ry)(ej • f,y )] \ri-rj\-^, where D is the 
magnitude of the dipole moments located at and Vj , 
is the direction of the dipole moment at r^, and Vij is the 
direction of the line connecting dipolcs, while jr^ — r^l is 
the distance between dipoles, with rj — {jxa,jya) being 
the lattice vector, and a being the lattice spacing. 

The bare Hamiltonian is now transformed into an ef- 
fective many-body Hamiltonian which includes the ef- 
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fects of screening, as the dipole interaction between two 
bare dipoles is renormalized (reduced) by the dielectric 
function erf(r,r'), leading to an effective screened dipole- 
dipole interaction V{r) = V^(r')e^^(r, r'). Although 
this effect is negligible at small filling factors, it may be 
substantial for filling factors u > f/3 as the dielectric 
function becomes sufficiently large to reduce the range of 
the interaction to a few neighbors. 

Since we are interested in superfluid phases, we choose 
to tune the experimental parameters to generate mostly 
attractive interactions. For the off-site dipolar inter- 
actions, we choose to have all dipoles aligned along 
the same direction {a, (p) of an external electric field, 
where a (0) is the polar (azimuthal) angle with re- 
spect to the z [x) axis. The interaction becomes Vij = 
Vx5(^.,j) {Vij = VyS(^ij-^) along the x (y) axis. Here, 
14 = D^{1 - 3 Sirica cos^(j>)/[a^edia)l Vy = D^{1 - 
3sin^ asin^ </i)/[a^ed(a)], while J^y) = 1 for nearest 
neighbors and zero otherwise. For the angles 4> = 
±7r/4, ±37r/4, the interactions are = Vy = V = 
D^{1 — Ssin^ a/2)/[a'^ed(a)], and become negative when 
the condition sin^ a > 2/3 is satisfied. For the on-site 
interaction, we choose to adjust the scattering length 
to produce U = -\U\ < 0. 

To study the physics discussed above, we use the effec- 
tive two-dimensional model hamiltonian 

H = -tY^ cj^Cj^ - \U\ rii^n.i - \V\ ^ n,„nj^, 

i {ij}<y<y' 

on a square lattice. In order to establish the quantum 
phases as a function of filling factor i/, we start by con- 
structing the partition hmction Z = j T)c)'Dce^ for the 
action 



S = I dT 
la 



(1) 

By symmetry, the singlet order parameters for super- 
fluidity correspond to local s-wave Asi, extended s-wave 
Age and d-wave pairing Upon integration of the 
fermionic degrees of freedom the action becomes 



q,a 



T T 



T 



where Ns is the total number of sites, a = si, se, d, the 
interactions Vsi = \U\, Vse = Vd = \V\, and the four- 
vector q = (w„,q). The inverse free fermion propagator 
matrix is 







iujn + 



3fe.fe' 



(3) 



with kinetic energy = fk — band dispersion 
Ek = — 2t [cos(fca;a) -I- cos(fcj,a)], chemical potential /x. 



four- vector k = (ia;„,k), and unit cell length a. The 
additional matrix appearing in Eq. ([2]) is 



V(fc,fc') = 







A„(fc - k') 
k') 



A„(k,k'), (4) 



where the Einstein summation over a is understood, and 
AQ(k, k') are the symmetry factors for the order param- 
eters, which in the limit of zero momentum pairing (k = 
k') become Asi(k, k) ~ I, Ase(k, k) = cos(fc2;a)+cos(fcj,a), 
Aci(k, k) ~ cos{kxa) — cos(fcj,a). 

In terms of the quasiparticle (7 = 2) or quasihole (7 = 
I) energies E'k,^ = (—)'•' \/^k + \^a^aO<-)\^ , where the 
symmetry function Ac((k) = AQ(k, k), the effective action 
becomes 



'\U\T 



\As 



\V\T 



(|A. 



\Ad\')+S2- 



IJ.Ns 
T 



(5) 



Here, the second term in the action is S2 = 
Sk 7 [1 + cxp(— i?k.7/r)]. Notice that there are three 
possible pure phases: local s-wavc {sl)\ extended s-wave 
(se) and c?-wave {d). In addition, there are several possi- 
ble binary mixed phases sZ ± se, si ±d and se ± d, which 
do not break time-reversal symmetry, and there are also 
those that do, such as si Use, sl±id and seHd. Lastly, 
several possible ternary mixed phases involving all three 
symmetries si, se and d may also exist. 

General Case: The order parameter equations can be 
obtained by minimization of the action with respective 
to each order parameter. By taking SS/SA'^ = 0, with 
a = si, se, d, we obtain 



Aa 



tanh(Sk,2/2r) 

N.- ^ 



2E^ 



k,2 



Aa(k), 



(6) 



with symmetry factors AQ,(k) — AQ,(k) [AQ'AQ,'(k)] , 
where repeated indices a' indicate summation. 

The number equation that fixes the chemical potential 
is obtained through the thermodynamic relation N = 
—dn/d^, where = —TlnZ is the thermodynamic po- 
tential. In the present approximation i7 = —TS, and the 
number equation reduces to 



-E 

Ns ^ 



1 - 
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k,2 



■tanh(£;k,2/2T) 



(7) 



where v = N/Ns is the filling factor. 

Using the amplitude-phase representation, we write 
the order parameters as Asi = |As;|e"^-"' for the local 
s-wave symmetry, A^e = |Ase|e*'^°'' for the extended s- 
wave symmetry, and A^ = |Arf|e"^'' for the c?-wave sym- 
metry. The filling factor dependence of the critical tem- 
perature Tc(i/) and the critical chemical potential /Xc(i^) 
can be obtained for pure si-, se-, and d-wavc symmetries 
by setting the order parameters A^; — 0, Ase — 0, and 
Ad = in Eqs. ^ and The solutions for Tc{v) are 
shown in Fig. [T]for two cases \U\/t = and \V\/t = 3, 
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and \U\lt = 2 and \V\/t = 3, where the corresponding 
supcrfluid phases are also indicated. The phase diagram 
is symmetric about v ^ 1, since the Helmholtz free en- 
ergy = fi + is invariant under the global particle- 
hole transformation ^ — > — and v ^ 2 ~ v. Notice 
that s-wave phases are favored at lower filling factors, 
while the d-wave phase is favored near half-filling, and 
this is directly correlated with the higher effective density 
of states in this vicinity. The time-reversal-symmetry- 
breaking phases occur at filling factors between the s- 
wave and d-wave phases. Generally, when \U\lt = 0, 
the only accessible phases are se-, d- and (se ± id)-wa.ve, 
however when \U\/t ^ the only accessible phases are 
{si + se)-, d and {si -I- se ± id)-wa.ve. 
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FIG. 1: Critical temperature Tc/t versus filling factor f at 
fixed interaction \U\/t = in (a) or \U\/t = 2 in (b) and 
\V\/t = 3. Notice the tetracritical point where the normal 
and all superconducting phases meet. 

Given that on-site interactions can be experimentally 
controlled, we focus our discussion at \U\/t = 0, which al- 
ready contains the essential physics of superfluid phases 
that spontaneously break time reversal symmetry and 
have a d-wave component. The Ginzburg-Landau theory 
near Tc is obtained by expanding the action of Eq. ([5]) 
in terms of the order parameters Age, and their 
complex conjugates. From the thermodynamic potential 
il = —TS, we can calculate the Helmholtz Free energy 
F = n + . The free energy per site F = F/Ng takes 
the simple form 



2b,d [1 + i cos(2(5</.)] |A,eP|Arf|2 + ^(i. - 1). 



(8) 



The cocfhcicnts a and b depend explicitly on the parame- 
ters of the model used and S(j) = (pd — '/'se- In the present 
case the possible phases se ± d are not accessible, and a 
tetracritical point exists where the normal and supercon- 
ducting phases with se, d and se ± id symmetries meet. 
In addition, the free energy depends only on cos(2(50) 
and does not distinguish between the phases se + id and 
se — id, which are thus degenerate. In the se ± id phases, 
timc-rcvcrsal symmetry is broken but not chirality. 

Harmonic Trap: The essential effect of an underlying 
harmonic trap Vh{r) — fcr^/2 is to allow for the emer- 
gence of non-uniform solutions. In particular, the har- 
monically confining potential allows for the emergence 



of all accessible phases se-, d- and (se ± zc?)-wave for 
\U\/t ~ and si + se, d- and {si -I- se ± id)-wave for 
\U\/t ^ 0. Within the local density approximation, we 
solve the order parameter Eq. ^ and number Eq. ([7]) 
with /i — > /i — V/i(r). In Fig. [21 we show the spatial profiles 
of the filling factor and order parameters as a function of 
dimensionless position rj = (i/e;j)^/^(r/a) measured from 
the center of the trap, where ~ ka^ /2. The parame- 
ters used arc \U\/t ~ 0, \V\/t = 3, T/t — 0.1, assuming 
that V = 1 (half-filling) at the center of the trap. Notice 
that as the filling factor decreases from the center of the 
trap to its edge, all accessible phases emerge: d-wave su- 
perfluid at the center of the trap, followed sequentially by 
regions of (seiid)- and se-wave superfluids, and the nor- 
mal state. Similarly, in the case of 7^ at low tem- 
peratures and assuming that v = 1 aX the center of the 
trap, the sequence of phases from the center of the trap is 
d-, {si -|- se ± id)-, {si + se)-wave superfluid followed by a 
normal region at the edge. The interesting qualitative as- 
pect here is the emergence of regions where time-reversal 
symmetry is spontaneously broken: (se±id) for \U\/t ~ 
and si -\- se zL id for \U\/t 7^ 0. This is very important 
in a very broad sense, because there are no confirmed 
examples in condensed matter physics of superfluid that 
spontaneously break time-reversal-symmetry. 
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FIG. 2: Spatially resolved filling factor u in (a) and super- 
fluid order parameters (a = se, d) in (b) as a function of 
77 = {t/eh)^^^{r/a), assuming that i/{0) — 1 (half-filling) at 
the center of the trap. The parameters used are \U\/t — 0, 
\V\/t = 3, and T/t = 0.1. 



Spontaneous Currents: In order to keep the discussion 
simple, we continue to focus on the case of \U\/t = 0, 
and discuss the spontaneous current flow in the shell cor- 
responding to the (se ± ic?)-wave superfluid. Consider 
for example that either the se -|- id phase or the se — id 
phase is realized in the example of Fig. [51 Given that ei- 
ther chiral phase spontaneously break time-reversal sym- 
metry, it is expected that within the boundaries of the 
se -f id {se — id) phase spontaneous currents circulate 
clockwise (counter-clockwise) near the outer boundary, 
and counter-clockwise (clockwise) near the inner bound- 
ary. To visualize the spontaneously generated currents, 
wc perform a long-wavelength expansion of the action in 
Eq. (|T]), which leads to the effective Free energy density 
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J'cS = J'di + J^nd + J^h + J^- The first term is 



-VA,, 



VA: 



: Cd,d 



VArf 



2m °2m 
the second term is non-diagonal in the indices se and d 



aA! 



2m 



2m 



■a,A<i + c.c. 



the third term is = 7seV/i(r)|As 
while the last term T is given in 
a current source term —idm — fti 



|2+7,y,(r)|A,|2, 
Eq. dHl). Adding 
and considering 
the phase difference 64> ~ 4'd — 4>se ~ we 
obtain within the se ± id phase the particle current 
density Ji = Ji^^ + ^i,|A|i in Cartesian representation 
{i = x,y). Here, Jj,^ = ^ [|AdpCd,d+ lA^epCse.^e] Sj^d 
is a phase-related contribution and Ji,\/\\ = 
[|^se|ci,se,d9j|Arf| - |Arf|cj,se,d9j|Ase|] is an 
amplitude-related component, where x = sm[5(j)) ~ ±1 
is the chirality of the seztid phase. In addition, the coeffi- 
cients Ct,se,d satisfy the relation Cx,se.d = -Cy.se.d = c^e.rf- 
Given the existence of the harmonic potential, we 
transform the currents to polar coordinates (?', 0), and 
require the radial current Jr to vanish ( = f • J = 0) 
at the boundaries between the se ± id and se occurring 
at r = Rge and at the boundaries between se ± id and d 
phases occurring at r = Rd- At these boundaries sponta- 
neous currents flow only within the se ± id phase limits, 
since these are the only phases that break spontaneously 
time-reversal symmetry. Under these conditions, a 
simple solution of the form (/id ~ xi'^/'^ ''r f{r)d] 
{cj)se = Xf{r)e) with f{r = Rse) = +1, /(r = Rd) = -1, 
and df{r)/dr\ii^^ = dJ{r)/dr\R^ = is pos- 
sible. The spontaneous currents at the inter- 
face boundaries are tangential, have the simple 
forms Jg(r = Rse) « {2x/m)\Ad\'^Cd,d/Rse and 
Je{r = Rd) ~ ~(2x/m)\Ase\'^Cse,se/Rd, and can be 
detected via Bragg spectroscopy as discussed next. 

Detection of time-reversal-symmetry-breaking: A de- 
tection scheme of spontaneous currents using velocity 
sensitive Bragg spectroscopy [13, [l^ is shown in Fig. |3] 
with right- (left-) going beam of frequency uj (w') and 
linear momentum k (k'). In Fig. [3^, circulating currents 
are shown at the boundaries of the region for (se-\-id) su- 
perfluidity, due to spontaneous breaking of time reversal 
symmetry at lower temperatures. The case of a normal 
region (higher temperatures), where no spontaneous cur- 
rents exist, is shown in Fig.ISjj for comparison. 

Conclusions: We discussed screened dipolar fermions 
in harmonically confined optical lattices modeled by an 
extended attractive Hubbard model, where both interac- 
tions and filling factors can be controlled. We established 
the superfluid phases in the singlet channel and indicated 
that accessible phases have not only pure s-wave or d- 
wave characters, but also mixed (s ± ic?)-wave character 



which breaks time reversal symmetry spontaneously. We 
calculated the spatially-dependent proflles of filling fac- 
tor and order-parameter for various supcrfiuid phases, 

y 






y 












' Jj ^ 






(b) 





FIG. 3: Schematic plots of a velocity sensitive Bragg spec- 
troscopy scheme to detect circulating currents of superfluid re- 
gions that break time reversal symmetry spontaneously. The 
case illustrated corresponds to \U\/t = with filling factor 
V = 1 at the center of the trap. 



and proposed a Bragg spectroscopy experiment to detect 
the time-reversal symmetry breaking phase, which con- 
tains spontaneously circulating supercurrents. 
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